Exogenous Brassinosteroid Enhances Zinc tolerance by activating the Phenylpropanoid Biosynthesis pathway in Citrullus lanatus L

ABSTRACT Zinc (Zn) is an important element in plants, but over-accumulation of Zn is harmful. The phytohormone brassinosteroids (BRs) play a key role in regulating plant growth, development, and response to stress. However, the role of BRs in watermelon (Citrullus lanatus L.) under Zn stress, one of the most important horticultural crops, remains largely unknown. In this study, we revealed that 24-epibrassinolide (EBR), a bioactive BR enhanced Zn tolerance in watermelon plants, which was related to the EBR-induced increase in the fresh weight, chlorophyll content, and net photosynthetic rate (Pn) and decrease in the content of hydrogen peroxide (H2O2), malondialdehyde (MDA), and Zn in watermelon leaves. Through RNA deep sequencing (RNA-seq), 350 different expressed genes (DEG) were found to be involved in the response to Zn stress after EBR treatment, including 175 up-regulated DEGs and 175 down-regulated DEGs. The up-regulated DEGs were significantly enriched in ‘phenylpropanoid biosynthesis’ pathway (map00940) using KEGG enrichment analysis. The gene expression levels of PAL, 4CL, CCR, and CCoAOMT, key genes involved in phenylpropanoid pathway, were significantly induced after EBR treatment. In addition, compared with Zn stress alone, EBR treatment significantly promoted the activities of PAL, 4CL, and POD by 30.90%, 20.69%, and 47.28%, respectively, and increased the content of total phenolic compounds, total flavonoids, and lignin by 23.02%, 40.37%, and 29.26%, respectively. The present research indicates that EBR plays an active role in strengthening Zn tolerance, thus providing new insights into the mechanism of BRs enhancing heavy metal tolerance.


Introduction
Heavy metals (HMs) pollution is a major environmental problem around the world. HMs are transported to the aboveground part through roots and accumulated in the edible part of crops, and enter human or animal bodies through the food chain, threatening human life and health 1 . Zinc (Zn), an essential biological element, participates in plant physiological metabolism as well as in the process of protein synthesis and transport 2 . However, high concentration of Zn will cause serious environmental pollution, and excessive accumulation of Zn in plants causes serious toxicity to plants 2,3 . Excessive Zn can cause the imbalance of nutrient and water absorption of plants, change the ultrastructure of chloroplasts to inhibit plant photosynthesis, destroy the membrane system and macromolecular substances of plant cells through accumulation of reactive oxygen species (ROS), and then seriously damage the growth, development, metabolism, and other plant physiological processes, causing plant death 3,4 . For example, high concentration of Zn in M. sativa produces a large amount of ROS, resulting in oxidative stress and cell damage, which seriously interferes M. sativa growth and quality 5 . Therefore, it has become a research hotspot in the field of environmental science to control Zn pollution, study the mechanism of plant tolerance to Zn and reduce the excessive absorption of Zn by plants.
Brassinolides (BRs) is a kind of natural plant hormone widely existing in plants. It has been identified as the sixth kind of plant hormone and plays a crucial role in plant growth and development, as well as biotic and abiotic stresses 6,7 . The regulation of exogenous BRs on HMs tolerance has attracted more and more attention 8 . Previous studies have shown that BRs can eliminate ROS produced by Zn stress by combining with membrane proteins, thereby alleviating the damage effect of Zn on plants. BRs can also affect the absorption of Zn by plants by stabilizing the electrical properties of cell membranes and enzyme activities, thus reducing Zn toxicity 9 . For example, Exogenous EBR spraying can improve eggplant seedlings tolerance to Zn stress through enhancement of antioxidant enzyme activity, osmotic substance accumulation, and hormone metabolism balance 10 . Exogenous EBR further increases the levels of antioxidant enzymes, thus effectively eliminating ROS, alleviating the toxicity of Zn to soybean seedlings 11 . Recently, EBR application to M. sativa under Zn stress could reduce Zn accumulation, promote the response of antioxidant defense system, and actively regulate the mechanism of heavy metal detoxification 5 .
Phenylpropane metabolic pathway is an important secondary metabolic pathway in plants, and its metabolites phenolics, flavonoids, and lignin are closely related to plant stress resistance 12,13 . Phenylpropane metabolism pathway plays an antioxidant role directly or indirectly in plant tolerance to heavy metal stress, and can improve plant absorption of heavy metal ions and stress tolerance 14 . In recent years, some studies have shown that BRs can improve plant stress resistance by regulating phenylpropane metabolism. Seed priming with EBR changes the growth and phenylpropanoid pathway of flax to water deficit 15 . Leaf application of EBR alleviates rice salt stress by up regulating secondary metabolites produced via the phenylalanine biosynthesis pathway 16 . EBR plays an active role in Colletotrichum fructicola resistance through the induction of lignin synthesis in Camellia sinensis 17 . Up to now, there are few reports on the research of heavy metal resistance mediated by BRs through the regulation of phenylpropane metabolism pathway, let alone Zn stress.
In the past decade, RNA-sequencing (RNA-seq) technique, a high-throughput method, has been used to investigate the global expression profiles 18 , which reveal the signal transduction pathways involved in Zn stress tolerance in some plants, such as wheat 19 , barley 20 , and Arabidopsis 21 . However, the specific gene expression profile of plant Zn tolerance induced by BR is still unclear. Watermelon (Citrullus lanatus) is one of the most important horticultural crops. However, so far, the possibility of using EBR to improve Zn tolerance of watermelon has not been studied, let alone the regulatory mechanism. Considering the increase of Zn content in horticultural soil and the stress improvement characteristics of BRs, the molecular mechanism of EBR regulating the response of watermelon to Zn stress were studied in the current study. 0.05 µM EBR could effectively alleviate Zn stress in watermelon. The gene expression pattern regulated by EBR under Zn stress was studied by RNA-seq and the different expressed genes (DEG) was significantly enriched in 'phenylpropanoid biosynthesis. In addition, the important enzymes and secondary metabolites involved in the phenylpropanoid biosynthesis pathway have also been proved to be related to EBR application. The results of this study will provide new insights into the mechanism of BR enhancing Zn tolerance in plants.

Stress treatment
The Zn (ZnSO 4 ·7 H 2 O) concentration was selected based on preliminary experiment, using 1.0, 2.5, 5.0, or 10.0 mM of Zn. It was found that plant growth was inhibited at 5 mM Zn but not completely inhibited (IC50) (data not shown). The EBR concentrations were selected using 0.025, 0.05, 0.10, 0.20, or 0.50 μM EBR and 0.05 μM EBR were selected.
When the plant reached the 4-leaf stage, healthy seedlings with consistent morphology were randomly divided into two groups: control and Zn treatment group. For EBR concentration screening, the following experimental design was carried out: plants were subjected to Zn and EBR-free nutrient solution (Control), or subjected to 5.0 mM Zn with different concentration of EBR (0, 0.025, 0.05, 0.10, 0.20, or 0.50 μM EBR). Seedlings were sprayed with EBR one day in advance and were watered with the nutrient solution containing Zn every 3 days. After 10 days of Zn/EBR application, the third leaves were frozen in liquid N 2 and kept at−80°C until analysis.

Photosynthetic pigments and net photosynthetic rate (Pn)
The 0.1 g fresh leaf samples were homogenized in 10 ml of 80% frozen acetone at 4°C at 10, 000 g centrifuge for 10 min. The absorbances of chlorophyll a and b in the supernatant at 663, 645, and 480 nm were recorded.
The Li-6400 portable photosynthetic instrument (LI-COR, Lincoln, Nebraska, USA) was used to measure Pn of the third true leaf of watermelon, and the instrument automatically recorded the net photosynthetic rate. The measured light intensity is 1,000 µmol m −2 s −1 , leaf temperature is 25°C, CO 2 concentration is 400 µmol mol −1 , and relative humidity is 75%.

Content measurement of Zn, H 2 O 2 , and MDA
For content measurement of Zn, the dried material was weighed, ground to a powder, and digested with a 1:3 mixture of HCl:HNO3. The digests were then dissolved in ultrapure water. Then, the digested sample were analyzed using a flameatomic absorption spectrometer (AAS, PerkinElmer) and the content was expressed as mg g −1 dry weight. H 2 O 2 and MDA content were measured according to the method of Wu et al. 22 .

RNA sequencing and differentially expressed genes (DEGs)
Total RNA from watermelon leaves of Zn+EBR and Zn treatments was isolated. After the RNA sample was qualified, the enrichment of mRNA and the construction of cDNA library were carried out. After the library inspection was qualified, the RNA-seq was carried out using an Illumina HiSeq TM 2000 platform. Three biological replicates were performed. Each sample got more than 6 Gb of data. After sequencing, the raw data was filtered to obtain high-quality clean data. The filtered clean reads were aligned to the reference genome (http://www.waterme londb.cn) using HISAT2 software 23 . Differentially expressed genes (DEGs) were analyzed using the DESeq2 package with | log2 Fold change| ≥ 1 and p-value<0.05. The fold change indicates the ratio of the expression amount between two groups.

GO and KEGG Enrichment
BLASTALL (v2.2.26) software was used to carry out annotations of GO, KEGG, COG, KOG, NR, Swissprot, Pfam databases to analyze the DEG functions. KEGG pathway enrichment was analyzed using the KEGG database (https:// www.kegg.jp/kegg/) 24 . Using R (v3.6.2) to make KEGG enrichment scatter diagram for the first 20 most significant pathways. Gene ontology (GO) enrichment of DEGs was annotated into GO database, enriching and analyzing the BP, MF, and CC functions, respectively, using Fisher exact test. GO and KEGG terms with corrected p-value<0.05 were considered significantly enriched.

Validation by RT-qPCR
Total RNA from watermelon leaves of Zn+EBR and Zn treatments was isolated and the SuperscriptIII first-strand synthesis system (Invitrogen, Shanghai, China) was used for cDNA synthesis according to the manufacturer's protocol. Use SYBR Premium Ex Taq (Takara, Dalian, China) to conduct RT-qPCR in real-time PCR system. The selected specific primers for DEG were designed according to the Watermelon Genome Database (http://www.water melondb.cn), as shown in Supplementary Table S1. β-Actin gene was used as a reference for normalizing the data.
Determination of the activities of phenylalanine ammonialyase (PAL), 4-coumaric ligase (4CL), and peroxidase (POD) Plant material was homogenized with 50 mM sodium borate buffer (pH 8.0) supplying 1 mM ethylenediaminetetraacetic acid disodium salt (EDTA) and 2% polyvinyl pyrrolidone (PVP), and then immediately centrifuged at 12,000 g and 4°C for 20 min. The supernatant fraction was frozen in liquid nitrogen and placed at−80°C for determining enzyme activities of PAL, 4CL, and POD according to the method of Shao et al. 25 . Mainly, the activities of PAL, 4CL, POD and CAD were determined by ELISA kit (JiangLai, Shanghai, China) following the protocol.

Total phenols, total flavonoids, and lignin content
According to the method of Alam et al. 26 , the phenol content was estimated by Folin -Ciocalteu reagent, and the phenol content was expressed as mg gallic acid equivalent (GAE) g −1 of sample. The flavonoid was determined by sodium nitrite aluminum nitrate colorimetry with sophorin as the standard sample. The flavonoid content was expressed as mg g −1 sample. The lignin content was determined by spectrophotometry according to the method of Xu et al. 27 . Acetyl lignin is produced after the phenol hydroxyl in lignin is acetylated. The product has a characteristic absorption peak at 280 nm. The lignin content was measured on the basis of the changing absorbance.

Statistical analysis
Each treatment has three biological replicates. The data is the average ± SD (Standard Deviation) of the replicates displayed by the vertical error bar. One-way analysis of variance (ANOVA) and the Least Significance Difference (LSD) test (significance level is 0.05, P value≤0.05) were used to analyze the differences. SPSS version 20.0 (IBM, USA) was used for statistical analysis.

Growth and Zn accumulation
To explore the effect of EBR on watermelon seedlings against Zn stress, a concentration range (0, 0.025, 0.05, 0.10, 0.20, or 0.50 μM EBR) was used. The growth of watermelon was obviously inhibited under Zn treatment alone (Figure 1a). Compared with the control, the shoot fresh weight decreased significantly. However, after pre-spraying EBR with different concentrations, the inhibition of watermelon by Zn was alleviated, which depended on the concentration effect of EBR. Compared with Zn alone, 0.025 and 0.05 significantly increased shoot fresh weight by 14.83% and 16.68%, respectively, while the growth of watermelon was inhibited by 13.82% under 0.50 μM EBR. The results showed that 0.025-0.05 EBR could alleviate Zn damage and promote seedling growth of watermelon.
EBR-induced watermelon tolerance to Zn stress is accompanied by a decrease in Zn accumulation (Figure 1b). Compared with the control, the leaves treated with Zn showed accumulation of large amount of Zn, but 0.05 μM EBR pretreatment significantly reduced Zn accumulation, presenting the best dose manner (Figure 1b).

Chlorophyll and photosynthesis
The Chl content and Pn in watermelon seedlings under Zn treatment alone decreased significantly compared with the control (Figure 1c, d). The application of low concentration EBR, however, promoted the increase of Chl content and Pn. Compared with Zn alone, under 0.05 μM EBR treatment, Chl content, and Pn increased by 32.94% and 38.84%, respectively, while 0.50 μM EBR treatment significantly decreased Chl content and Pn. It indicates that exogenous low concentration EBR is helpful to increase the Chl content of watermelon seedlings under Zn stress, improve the net photosynthetic rate, and thus enhance the photosynthetic capacity.

Oxidative stress
Compared with the control, Zn alone stress showed an increase in the levels of H 2 O 2 and MDA in leaves ( Figure 2). Under Zn stress, spraying exogenous EBR increased watermelon tolerance to Zn in a dose-dependent manner, and 0.05 μM was an optimum concentration for spraying watermelon seedlings to decrease oxidative stress. 0.05 μM EBR could significantly reduce the content of H 2 O 2 and MDA in leaves by 36.57% and 32.79%, respectively, indicating that EBR might improve Zn stress tolerance by alleviating oxidative stress damage. Based on the above growth and physiological data, 0.05 μM EBR as a beneficial dose was selected for the following analysis.

Differential expression genes (DEGs)
Six cDNA libraries from two treatments (Zn; Zn +0.05 μM EBR) were sequenced using the Illumina HiSeq TM 2000 platform. The raw data contained total average 7,134,478,450 baseSum. After the low-quality sequences were removed, total average 24,448,262 readSum was obtained. The average of Q20, Q30 values, and GC content was 95.9%, 91.6%, and 45.5%, respectively (Table S2). These reads were, thus, mapped to the watermelon reference genome (http://www.waterme londb.cn). The DEGs analysis was conducted with the screening threshold was: | log2FC| ≥1 and q-adjust<0.05. We identified 350 DEGs between Zn+EBR and Zn treatments (Figure 3a; Table S3) including 175 up-regulated and 175 down-regulated DEGs. Hierarchical clustering was used to observe gene expression patterns with log10 FPKMs in both groups (Figure 3b).

Functional enrichment analysis of DEGs
GO enrichment analysis was used to evaluate the function of DEG, revealing the BP, MF, and CC categories for the 350 DEGs (Table S4) To check the DEG-related pathways, they were searched using the KEGG pathway database. KEGG enrichment indicated that the up-regulated DEGs were significantly associated with 'phenylpropanoid biosynthesis (map00940)', 'protein processing in endoplasmic reticulum (map04141)', and 'biosynthesis of secondary metabolites -unclassified (map00999)' (Figure 4a; Table S5). For down-regulated DEGs, no pathway was significantly enriched (Figure 4b).

Critical functional DEGs regulated by EBR under Zn stress
DEGs related to oxidative stress. The biological process of upregulated DEGs was significantly enriched to 'response to oxidative stress, such as common antioxidant enzyme genes peroxidase (ClG42_01g0049000, ClG42_10g0201300, ClG42_02g0242500, ClG42_09g0159900) and polyphenol oxidase (ClG42_03g0058500) ( Table S4), suggesting that EBRalleviated Zn stress was related to EBR-induced antioxidant protection to alleviate oxidative damage, which was consistent with the physiological data (Figures 1, 2).
DEGs related to cell wall. The cellular component (CC) of upregulated DEGs was significantly enriched in extracellular region and cell wall (Table S4), such as ClG42_05g0027500 encoding pectinesterase 2 and ClG42_10g0192200, ClG42_09g0179200, ClG42_08g0082000, and ClG42_07g0166400 encoding extensinlike protein enriched in external encapsulating structure and cell wall organization. Besides, ClG42_01g0049000 is a lignin- forming anionic peroxidase. PAL, 4CL, CCR, and COMT are involved in lignin biosynthetic process.
DEGs related to protein processing in the endoplasmic reticulum (ER). There were 12 DEGs involved in protein processing in the ER ( Figure S2; Table S5). Among them, the 11 differentially expressed HSPs and a E3 ubiquitin-protein ligase RMA1 genes were up-regulated. The DEGs were all enriched in the degradation of ubiquitin -proteasome-dependent process of ER-associated degradation (ERAD). The results indicate that ERAD is involved in EBR mediated Zn stress tolerance.

EBR activates phenylpropanoid biosynthesis pathway
Due to the significant enrichment in the phenylpropanoid biosynthesis pathway (Figure 4c; Figure S1), here, the activities of three key enzymes including PAL, 4CL, and POD were analyzed. The presence of EBR significantly increased the activities of PAL, 4CL, and POD of watermelon leaves by 30.90%, 20.69%, and 47.28%, respectively, under Zn stress compared to that in the absence of EBR ( Figure 5). Therefore, we further measured the content of these metabolites. As shown in Figure 5, compared with Zn alone, EBR promoted the accumulation of these substances in watermelon leaves. Total phenolic compounds, total flavonoids, and lignin in Zn+EBR treated leaves were averagely 23.02%, 40.37%, and 29.26% more than those in Zn treated leaves, respectively.

Discussion
When a large amount of Zn is enriched in the soil, it will adversely affect the growth of plants, then affect the quality and yield 1 . Here, Zn stress treatment significantly inhibited the growth of watermelon, and also significantly inhibited the accumulation of Zn, chlorophyll content, and photosynthetic rate of watermelon (Figure 1), which was consistent with the previous studies that plant growth was inhibited under Zn stress 4,5,[28][29][30] . To reduce the serious impact of heavy metal stress (including Zn) on plants, many methods have been adopted, including foliar application of different chemicals/ phytohormones 4,5,29,[31][32][33][34][35] . BRs act as multidimensional regulators of plant responses under different environmental stresses, and epibrassinolide (EBR) is a synthetic brassinolide analog, which has been widely used in production 7 . In this study, the addition of EBR could significantly promote the increase of the biomass of watermelon seedlings under Zn stress (Figure 1), indicating that exogenous EBR treatment can alleviate the watermelon growth inhibition caused by Zn stress, and 0.05 μM EBR treatment is the best. Similar studies have been reported in other different plants [28][29][30] . When plants are damaged by heavy metals, on the one hand, EBR can increase the activation of H± ATPase and activate cell wall lysis enzyme, thus promoting cell growth through cell division and elongation; on the other hand, EBR can increase photosynthesis, which may be the main reason 36 . Photosynthetic pigments are involved in the absorption, transmission, and transformation of light energy in plants, and are one of the indicators of photosynthetic capacity. When plants are in a stress environment for a long time, there will be a decrease in chlorophyll content, photosynthesis rate and other characteristics 37 . In this experiment, chlorophyll content and Pn in watermelon leaves decreased under Zn stress (Figure 1), indicating that Zn-induced decline of biomass is directly related to the decline of net photosynthetic rate of leaves. Low concentration EBR application could significantly increase the chlorophyll content and Pn in watermelon leaves, suggesting that EBR could promote the synthesis of chlorophyll or inhibit chlorophyllase activity in watermelon leaves, thus promoting the absorption and utilization of light energy in watermelon leaves 38 . Besides, EBR-induced increase in the levels of chlorophyll under stress was related to the decrease of ROS accumulation, thus reducing oxidative damage to thylakoid membrane structure and function 34 . The possible antioxidant defense mechanism for the beneficial effect of EBR application was discussed in the following sections.
Cell membrane is an important structure for plant cells to isolate protoplasts from the external environment, which can control the material exchange and signal transmission between cells and the external environments 39 . When plants are under stress, plant cells will produce excessive ROS, which will cause changes in plant membrane and damage to membrane structure. The plasma membrane and the inner membrane system of various organelles will expand or damage 39,40 . The accumulated active oxygen radicals will cause the unsaturated bonds in membrane fatty acids to be oxidized to form MDA. The level of MDA is an important indicator to reflect the strength of cells under optimal stress and the degree of plasma membrane damage 30 . Similar to other environmental stresses, excessive Zn leads to the formation of excessive ROS in plants cells resulting in plant cell oxidative damage and membrane lipid peroxidation 3,41 . This study showed that the H 2 O 2 and MDA content of watermelon seedlings were significantly increased after being stressed by Zn (Figure 2), indicating that the cell plasma membrane had peroxidation, which disturbed the normal physiological function of the plasma membrane. However, after exogenous EBR treatment, the content of H 2 O 2 and MDA decreased significantly, indicating that EBR had a positive role in mitigating oxidative stress damage caused by Zn stress. The biological process of up regulating DEG significantly enriched the response to oxidative stress (Table S4), which further demonstrated the mechanism of EBR mediated oxidative stress response to alleviate Zn damage. The main way for plant systems to overcome the oxidative stress is to use primary metabolites and/or secondary metabolites 42 . On the one hand, ROS neutralization is achieved through the complex antioxidant system developed by plant cells. Stress tolerance of major metabolites is mainly through classical antioxidant enzyme system and non-enzyme system 42 .
Here, under Zn stress exogenous EBR treatment could significantly increase the expression level of three peroxidase and one polyphenol oxidase genes (Table S4) and have a high activity of POD ( Figure 5). Peroxidases as 'stress enzymes' are heme containing proteins that use H 2 O 2 to oxidize various organic and inorganic substrates. Peroxidase activity is a sensitive indicator of heavy metal stress 43 . Polyphenol oxidase (PPO) acts as a quencher of chloroplast photooxidation. In addition, due to the vacuole position of PPO substrate, PPO may synthesize specific metabolites and react to environmental factors 44 . It is also believed that the polymerization of polyphenols by peroxidase increases after heavy metal absorption and detoxification, which is the reason for heavy metal binding in the skin gland of Nymphaea epidermal 45 . Alam et al. 26 found that providing EBR to soybean plants under NaCl stress led to the enhancement of peroxidase activities coupled with the reduction in H 2 O 2 and MDA content. These results indicated that exogenous EBR with appropriate concentration can stimulate and induce the antioxidant enzyme protection system, at least including the gene expression and activity of POD, accelerate the clearance of ROS, reduce ROS accumulation in cells, and reduce membrane lipid peroxidation damage.
On the other hand, the secondary defense also plays important roles in scavenging ROS by using a variety of secondary metabolites synthesized in plants like lignin, tannin, flavonoids, phenols, etc. 42 . These compounds are important metabolic intermediates of phenylpropane metabolic pathway. In the present study, KEGG enrichment revealed that the upregulated DEGs were significantly associated with 'phenylpropanoid biosynthesis (map00940)', and 'biosynthesis of secondary metabolites-unclassified (map00999)' (Figure 4a; Table S5), suggesting that EBR may help to reduce the Zn toxicity in watermelon through this pathway. The genes encoding enzymes in phenylpropanoid biosynthesis including PAL, 4CL, CCR, CALDH, and CAD, which are involved in the synthesis of flavonoids, isoflavones, and lignin, may be involved in Cd tolerance 46 . The phenylpropanoid biosynthesis pathway starts from phenylalanine, then PAL catalyzes the conversion of phenylalanine to trans-cinnamic acid. This pathway can form other phenolic compounds like flavonoids, coumarins, lignans, and lignin. They have antioxidant characteristics and can play an important role in ROS removal and the formation of metal complexes to protect plants from abiotic stress 12 . Through a series of reactions catalyzed by 4CL, C4H, HCT, and CCoAOMT, coumarinyl-CoA is converted into p-coumaroyl-CoA, caffeoyl-CoA, and feruloyl-CoA, which are the core modules of the phenylpropanoid pathway 13 . Peroxidase catalyzes the oxidation of phenylpropanoids to phenoxy, and the subsequent nonenzymatic coupling controls the pattern and extent of polymerization 47 . Under lead stress, phenylpropanoid biosynthesis is most significantly enriched in T. orientalis and the expression of DEGs such as PAL, C4H, 4CL and CCoAOM is up-regulated 48 . Bacillus altitudinis WR10 regulates the phenylpropanoid biosynthesis related genes, which may improve phenolic acids accumulation, thus protecting plant cells from copper toxicity 49 . In this study, EBR treatment up-regulated the expression levels of some key genes which were involved in phenylpropanoid biosynthesis under Zn stress, including PAL, 4CL, CCR, and CCoAOMT (Figure 4c; Figure S1), and the activities of three key enzymes like PAL, 4CL, and POD ( Figure 5), playing important roles in lignin, flavonoids, and phenolic compounds. The changes of content of lignin, flavonoids, and phenolics are consistent with those of these genes and enzymes ( Figure 5). All these results indicated that phenylpropanoid biosynthesis is an important module involved in EBR-induced Zn tolerance.
It is speculated that there are two mechanisms for phenylpropanoid pathway to participate in EBR-mediated Zn tolerance: one is that EBR further activates phenylpropyl biosynthesis pathway, leading to the generation of various phenolic compounds, which may eliminate harmful ROS; the other is to block Zn transport in plants. The rich metabolites produced by phenylpropane biosynthesis pathway coupled with high antioxidant capacity give mulberry higher salt tolerance 50 . Flavonoids are the largest class of special phenolic compounds synthesized through phenylpropane pathway, which can improve the tolerance of different plants by minimizing oxidative damage 16 . EBR supplementation further enhanced the flavonoid content in Glycine max 26 and Camellia sinensis 51 , and Vitis vinifera 52 and Oryza sativa 16 . Phenolic compounds are powerful heavy metal chelators. The preferential activation of genes involved in the early phenylpropanoid pathway may also be consistent with the increase of single molecule and final lignin production. EBR activates map00940-related pathways, increases lignin content in the secondary cell wall, which is closely associated with Zn absorption and transport, thus making EBR-treated plants more tolerant. At the same time, this is consistent with the enrichment of DEGs in the extracellular region and cell wall (Table S4).

Conclusion
The present study showed that exogenous EBR enhanced Zn tolerance in watermelon seedlings with higher levels of chlorophyll and Pn and lower accumulation of H 2 O 2 and MDA in watermelon leaves. The comparative transcriptome data and physiological analysis revealed that the increased H 2 O 2 scavenging and activated phenylpropanoid biosynthesis pathway contributing to the mitigation of Zn toxicity by EBR. As far as we know, this is the first study to clarify the mechanism of BR reducing Zn stress through phenylpropanoid biosynthesis pathway. In conclusion, these results enhance our understanding of the mechanism of exogenous EBR enhancing Zn tolerance and the heavy metal responsive genes identified after EBR application, providing a target for future molecular breeding. 53